FIRST EVIDENCE OF NEW PHYSICS IN b <-> s TRANSITIONS 
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We combine all the available experimental information on B 3 mixing, including the very recent 
tagged analyses of B s — * J /^<f> by the CDF and D0 collaborations. We find that the phase of 
the B 3 mixing amplitude deviates more than 3a from the Standard Model prediction. While no 
single measurement has a 3cr significance yet, all the constraints show a remarkable agreement with 
the combined result. This is a first evidence of physics beyond the Standard Model. This result 
disfavours New Physics models with Minimal Flavour Violation with the same significance. 



In the Standard Model (SM), all flavour and CP vio- 
lating phenomena in weak decays are described in terms 
of quark masses and the four independent parameters 
in the Cabibbo-Kobayashi-Maskawa (CKM) matrix 
In particular, there is only one source of CP violation, 
which is connected to the area of the Unitarity Trian- 
gle (UT). A peculiar prediction of the SM, due to the 
hierarchy among CKM matrix elements, is that CP vio- 
lation in B s mixing should be tiny. This property is also 
valid in models of Minimal Flavour Violation (MFV) Q , 
where flavour and CP violation are still governed by the 
CKM matrix. Therefore, the experimental observation of 
sizable CP violation in B s mixing is a clear (and clean) 
signal of New Physics (NP) and a violation of the MFV 
paradigm. In the past decade, B factories have collected 
an impressive amount of data on Bd flavour- and CP- 
violating processes. The CKM paradigm has passed un- 
scathed all the tests performed at the B factories down 
to an accuracy just below 10% [H, 0|. This has been 
often considered as an indication pointing to the MFV 
hypothesis, which has received considerable attention in 
recent years. The only possible hint of non-MFV NP is 
found in the penguin-dominated b — > s non-leptonic de- 
cays. Indeed, in the SM, the S q q S coefficient of the time- 
dependent CP asymmetry in these channels is equal to 
the S cEs measured with b — > ccs decays, up to hadronic 
uncertainties related to subleading terms in the decay 
amplitudes. Present data show a systematic, although 
not statistically significant, downward shift of S q q S with 
respect to S cSs J>\ , while hadronic models predict a shift 
in the opposite direction in many cases [g, [7| . 

From the theoretical point of view, the hierarchical 



structure of quark masses and mixing angles of the SM 
calls for an explanation in terms of flavour symmetries 
or of other dynamical mechanisms, such as, for example, 
fermion localization in models with extra dimensions. All 
such explanations depart from the MFV paradigm, and 
generically cause deviations from the SM in flavour vio- 
lating processes. Models with localized fermions Q, and 
more generally models of Next-to-Minimal Flavour Vio- 
lation [9| , tend to produce too large effects in e K 10, 11}. 
On the contrary, flavour models based on nonabelian 
flavour symmetries, such as U (2) or SU (3), typically sup- 
press NP contributions to s <-> d and possibly also to 
b «-> d transitions, but easily produce large NP contribu- 
tions to b <-> s processes. This is due to the large flavour 
symmetry breaking caused by the top quark Yukawa cou- 
pling. Thus, if (nonabelian) flavour symmetry models are 
relevant for the solution of the SM flavour problem, one 
expects on general grounds NP contributions to b <-> s 
transitions. On the other hand, in the context of Grand 
Unified Theories (GUTs), there is a connection between 
leptonic and hadronic flavour violation. In particular, in 
a broad class of GUTs, the large mixing angle observed 
in neutrino oscillations corresponds to large NP contri- 
butions to b <-> s transitions [12| ■ 

In this Letter, we show that present data give evidence 
of a B s mixing phase much larger than expected in the 
SM, with a significance of more than 3a. This result is 
obtained by combining all available experimental infor- 
mation with the method used by our collaboration for 
UT analyses and described in Ref. [l3j |. 

We perform a model-independent analysis of NP con- 
tributions to B s mixing using the following parametriza- 
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where H^g 1 is the effective Hamiltonian generated 
by both SM and NP, while only contains SM 

contributions. The angle j3 s is defined as f3 s = 
axg(-(V ts V t l)/(V cs V* b )) and it equals 0.018 ± 0.001 in 
the SM. 1 
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0.19 ± 0.07 
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1.52 ± 0.06 
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TABLE I: Input parameters used in the analysis. We also 
show the correlation coefficients Cs of the measurements of 
A s , AF S and tr from ref. 1211 1 . 
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We use the followin g ex perimental input: the CDF 
measurement of Am s I15l , the semileptonic asymme- 
try in B s decays Aqt (l6j . the dimuon charge asymme- 
try from D0 17J and CDF [l8j, the measurement 
of the B s lifetime from flavour-specific final states [l9j |. 
the two-dimensional likelihood ratio for Ar s and 4> s = 
2(/3 s —(pB a ) from the time-dependent tagged angular anal- 
ysis of B s — > J I ip(f> decays by CDF [2(| and the correlated 
constraints on r s , Ar s and <j) s from the same analysis 
performed by D0 [2lj. For the latter, since the com- 
plete likelihood is not available yet, we start from the 
results of the 7-variable fit in the free-^., case from Ta- 
ble I of ref. 
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We implement the 7x7 correlation 
matrix and integrate over the strong phases and decay 
amplitudes to obtain the reduced 3x3 correlation ma- 
trix used in our analysis. In the D0 analysis, the twofold 
ambiguity inherent in the measurement (<j) s — > 7r — cf> s , 
AT S — » — Ar s , cos — > — cos £1,3) for arbitrary strong 
phases was removed using a value for cos <5i ; 2 derived from 
the BaBar analysis of B d — > J/^K* using SU(3). How- 
ever, the strong phases in Bd — > J/^K* and B s — > 
cannot be exactly related in the SU(3) limit due to the 



1 We are using the usual CKM phase convention in which VtsV^ 
is real to a very good approximation. 
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FIG. 1: From left to right and from top to bottom, 68% 
(dark) and 95% (light) probability regions in the Ab„ - Cb s , 
A^ p /A s a M - AT planes and p.d.f for C Bs , A B „ A^ P /A S 3 M , 
- Re A^/A™, Im y£ P Mf M . 



/NP 



singlet component of <j>. Although the sign of cos^i.2 ob- 
tained using SU(3) is consistent with the factorization es- 
timate, to be conservative we reintroduce the ambiguity 
in the D0 measurement. To this end, we take the errors 
quoted by D0 as Gaussian and duplicate the likelihood 
at the point obtained by applying the discrete ambigu- 
ity. Indeed, looking at Fig. 2 of ref. 21 1, this seems 
a reasonable procedure. Hopefully D0 will present re- 
sults without assumptions on the strong phases in the 
future, allowing for a more straightforward combination. 
Finally, for the CKM parameters we perform the UT 
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analysis in the presence of arbitrary NP as described in 
ref. [nj, obtaining p = 0.140 ± 0.046, fj = 0.384 ± 0.035 
and sin2/3 s = 0.0409 ± 0.0038. The new input param- 
eters used in our analysis are summarized in Table [H 
all the others are given in Ref. The relevant NLO 

formulae for Ar s and for the semilcptonic asymmetries 
in the presence of NP have been already discussed in 
refs. 
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FIG. 2: P.d.f. for 0s 3 without the tagged analysis of B s — > 
J/^<f) (top left), including only the CDF analysis (top right), 
including only the D0 Gaussian analysis (bottom left), in- 
cluding only the D0 likelihood profiles (bottom right). We 
show 68% (dark) and 95% (light) probability regions. 



TABLE II: Fit results for NP parameters, semileptonic asym- 
metries and width differences. Whenever present, we list the 
two solutions due to the ambiguity of the measurements. The 
first line corresponds to the one closer to the SM. 

The results of our analysis are summarized in Table 
Ull We see that the phase 4>b s deviates from zero at 
3.7a\ We comment below on the stability of this signif- 
icance. In Fig. Q] we present the two-dimensional 68% 
and 95% probability regions for the NP parameters Cb. 
and 4>b 3 i the corresponding regions for the parameters 
A™ P /A^ M and </>^ p , and the one-dimensional distribu- 
tions for NP parameters. Notice that the ambiguity of 
the tagged analysis of B s — > J/^<f> is slightly broken by 
the presence of the CKM-sublcading terms in the expres- 
sion of T12/M12 (see for example eq. (5) of ref. t 23]). The 
solution around 4>b s ~ —20° corresponds to ~ —50° 
and A^ p /Ag ~ 75%. The second solution is much more 
distant from the SM and it requires a dominant NP con- 
tribution (A^ P /A^ M ~ 190%). In this case the NP phase 
is thus very well determined. The strong phase ambigu- 
ity affects the sign of cos <j) s and thus Re A^ p /A^ M , while 
Im Af F /A^ M 0.74 in any case. 

Before concluding, we comment on our treatment of 
the D0 result for the tagged analysis and on the sta- 
bility of the NP fit. Clearly, the procedure to reintro- 
duce the strong phase ambiguity in the D0 result and 
to combine it with CDF is not unique given the avail- 
able information. In particular, the Gaussian assumption 
can be questioned, given the likelihood profiles shown in 



Ref. [2l[. Thus, we have tested the significance of the NP 
signal against different modeling of the probability den- 
sity function (p.d.f.). First, we have used the 90% C.L. 
range for <j) s = [—0.06, 1.20]° given by D0 to estimate the 
standard deviation, obtaining cf> s — (0.57 ± 0.38)° as in- 
put for our Gaussian analysis. This is conservative since 
the likelihood has a visibly larger half-width on the side 
opposite to the SM expectation (see Fig. 2 of Ref. [2lj). 
Second, we have implemented the likelihood profiles for 
(j) s and AT S given by D0, discarding the correlations 
but restoring the strong phase ambiguity. The likelihood 
profiles include the second minimum corresponding to 
4> s — > 4> s + 7T, Ar — > — Ar, which is disfavoured by the 
oscillating terms present in the tagged analysis and is dis- 
carded in our Gaussian analysis. Also this approach is 
conservative since each one-dimensional profile likelihood 
is minimized with respect to the other variables relevant 
for our analysis. It is remarkable that both methods give 
a deviation of <pB a from zero of 3 u (the 3 u ranges for <pB a 
are [-88, -48]°U[-41, 0]° and [-88, 0]° for the two meth- 
ods respectively). We conclude that the combined anal- 
ysis gives a stable evidence for NP, although the precise 
number of standard deviations depends on the procedure 
followed to combine presently available data. 

To illustrate the impact of the experimental con- 
straints, we show in Fig. [2] the p.d.f. for (f>B B obtained 
without the tagged analysis of B s — > J/^><j> or includ- 
ing only CDF or D0 results. Including only the CDF 
tagged analysis, we obtain 4>b s < at 97.7% probability 
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(2.3a). For D0, we show results obtained with the Gaus- 
sian and likelihood profile treatment of the errors. In the 
Gaussian case, the D0 tagged analysis gives 4>b s < 
at 98.0% probability (2.3c), while using the likelihood 
profiles <p Bs < at 92.8% probability (1.8er). Finally 
it is remarkable that the different constraints in Fig. [5] 
are all consistent among themselves and with the com- 
bined result. We notice, however, that the top-left plot is 
dominated by the measurement of Ag£ while A| L favours 
positive 4>b 3 , although with a very low significance. For 
completeness, in Table [TTI we also quote the fit results for 
Af L , Ag£ and for Ar s /r s . 

In this Letter we have presented the combination of all 
available constraints on the B s mixing amplitude lead- 
ing to a first evidence of NP contributions to the CP- 
violating phase. With the procedure we followed to com- 
bine the available data, we obtain an evidence for NP at 
more than 3c To put this conclusion on firmer grounds, 
it would be advisable to combine the likelihoods of the 
tagged B s — > J/^cf) angular analyses obtained without 
theoretical assumptions. This should be feasible in the 
near future. We are eager to see updated measurements 
using larger data sets from both the Tevatron experi- 
ments in order to strengthen the present evidence, wait- 
ing for the advent of LHCb for a high-precision measure- 
ment of the NP phase. 

It is remarkable that to explain the result obtained for 
(f> s , new sources of CP violation beyond the CKM phase 
are required, strongly disfavouring the MFV hypothesis. 
These new phases will in general produce correlated ef- 
fects in AB = 2 processes and in b — ► s decays. These 
correlations cannot be studied in a model-independent 
way, but it will be interesting to analyse them in specific 
extensions of the SM. In this respect, improving the re- 
sults on CP violation in b — > s penguins at present and 
future experimental facilities is of the utmost importance. 



Note added 

During the review procedure of this Letter, results 
based on new data were presented by the Tevatron ex- 
periments, as well as a combination of Tevatron results 
on the tagged angular analysis of B s — > J/tptp. How- 
ever these updates are all unpublished. Furthermore, 
the likelihoods required by our analysis are not publicly 
available except for the new D0 analysis with no as- 
sumption on strong phases 



24|. For the sake of com- 
pleteness, we quote (j> Bs = (-19 ± 8)° U (-69 ± 7)° 
([-36,-5]° U [-83,-54]° at 95% probability), obtained 
using this new likelihood for the D0 tagged angular anal- 
ysis of B s — > J/ipcf). Clearly, we no longer need to ma- 
nipulate the D0 to remove the strong phase assumption 
and to account for the non-Gaussian shape as described 
above. Remarkably, this updated result is well compati- 
ble with the results of this Letter, confirming a deviation 



from the SM at the level of ~ 3cr (99.6% probability). 
More recent experimental results seem to confirm the ef- 
fect discussed in this Letter. We will include them in 
future analyses as soon as they become available. 
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this analysis and for improving it with several valuable 
suggestions. We also thank G. Giurgiu, G. Punzi and 
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perimental results. We acknowledge partial support from 
RTN European contracts MRTN-CT-2006-035482 "FLA- 
VIAnet" and MRTN-CT-2006-035505 "Heptools". M.C. 
is associated to the Dipartimento di Fisica, Universita di 
Roma Tre. E.F. and L.S. are associated to the Diparti- 
mento di Fisica, Universita di Roma "La Sapienza" . 
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